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The combination of General Relativity and the Standard Model of particle physics 
disagree with numerous observations on scales from our Solar System up. In the 
canoncial concordance model of ACDM cosmology, many of these contradictions 
between theory and data are removed or alleviated by the introduction of three 
completely independent new components of stress-energy - the inflaton, dark mat- 
ter, and dark energy. Each of these in its turn is meant to have (or to currently) 
dominate the dynamics of the universe. There is, until now, no non-gravitational 
evidence for any of these dark sectors; nor is there evidence (though there may be 
motivation) for the required extension of the standard model. An alternative is to 
imagine that it is General Relativity that must be modified to account for some or 
all of these disagreements. Certain coincidences of scale even suggest that one might 
expect not to have to make independent modifications of the theory to replace each 
of the three dark sectors. Because they must address the most different types of 
data, attempts to replace dark matter with modified gravity are the most con- 
troversial. A phenomenological model (or family of models), Modified Newtonian 
Dynamics, has, over the last few years seen several covariant realizations. We dis- 
cuss a number of challenges that any model that seeks to replace dark matter with 
modified gravity must face: the loss of Birkhoff's Theorem, and the calculational 
simplifications it implies; the failure to explain clusters, whether static or interact- 
ing, and the consequent need to introduce dark matter of some form, whether hot 
dark matter neutrinos, or dark fields that arise in new sectors of the modified grav- 
ity theory; the intrusion of cosmological expansion into the modified force law, that 
arises precisely because of the coincidence in scale between the centripetal acceler- 
ation at which Newtonian gravity fails in galaxies, and the cosmic acceleration. We 
conclude with the observation that, although modified gravity may indeed manage 
to replace dark matter, it is likely to do so by becoming or at least incorporating, 
a dark matter theory itself. 
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1. Introduction 

The Standard Model of particle physics, augmented by the generally benign addi- 
tion of a mechanism for generating neutrino flavor mixing, successfully describes 
a remarkable range of short distance phenomena. From the highest energy colli- 
sions at particle colliders (and in cosmic rays) to the internal structure of stars, 
observations match this model's predictions to the precision with which physicists 
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have been able to calculate them. Meanwhile from the scale of the Earth's orbit 
around the sun down to just a few microns, General Relativity accurately describes 
all known gravitational phenomena. 

Much has been made of our inability to reconcile quantum field theory - the 
mathematical and conceptual basis of the Standard Model - with general relativity 
on ultra-short scales. Addressing this failure is the primary motivation behind string 
theory, loop quantum gravity, and other approaches to quantum gravity. However, 
while the ultra-violet failings of the standard model and GR are as yet theoretical 
- no experiment currently probes the quantum gravity regime - the infrared fail- 
ings of this pair are real. At all scales above that of the solar system, the 
combination of General Relativity and the Standard Model contradict 
observations. 

Viewing the past century of astronomical and cosmological data as a massive 
assault on our canonical theories is a non-traditional perspective, but one that 
underlies the whole program of modified gravity. It is therefore appropriate to 
recall what a standard model general relativistic cosmology would look like and 
review the ways in which observations contradict its predictions. This we shall do 
in section one, where we shall also review the standard "dark side" solutions to 
these problems and see how it leads to a concordance cosmology that would have 
made William of Ockham shudder. 

In section two, we will present the possibility of modified gravity as a solution 
to the failure of general relativistic standard model cosmology, focusing on frame- 
works where it offers a comprehensive solution to both the missing mass in cosmic 
structures (aka the missing or dark matter problem) and the observation of cosmic 
acceleration (aka the dark energy problem). 

In the third section, we will discuss the challenges these comprehensive modified 
gravity theories face, and why they seem doomed to give Father William even less 
comfort than their dark competitors. 

2. The many failures of Standard Model General Relativistic 

Cosmology 

The standard model contains a very limited number of stable ingredients from 
which to assemble the universe we observe. Protons and neutrons (and their anti- 
particles) are the only stable baryons. They can be combined into just 255 known 
stable nuclei, ranging from hydrogen to lead-208, however, the universe contains 
only trace amounts of anything other than hydrogen ( 1 H in particular) and Helium- 
4. (Deuterium ( 2 H) has an abundance of 0.01% relative to 1 H, and Helium-3 an 
abundance of a few 10~ 5 . Among higher atomic number elements, oxygen has the 
greatest abundance, but is still less than 0.1%.) This paucity of heavy elements 
should come as something of a surprise since many nuclei are more stable than 
hydrogen or helium. It results from a variety of accidents of nuclear physics and 
cosmology. These include the combination of the small binding energy of deuterium 
and the large number of photons per baryon in the universe which result in the 
so-called deuterium bottleneck - the delay in the synthesis of elements heavier 
than hydrogen until the universe has cooled to a temperature where characteristic 
energies are far far below the characteristic activation energies of nuclear fusion 
processes. The absence of any stable isotope of atomic mass number eight, further 
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inhibits processing of anything beyond lithium and beryllium. All heavier isotopes 
are produced almost exclusively much later - in stars and in interstellar cosmic 
ray interactions. Meanwhile, the only stable negatively charged particle that can 
neutralize the charge of a nucleus to form an atom is an electron. 

From the point of view of a cosmologist, this parsimony is fortunate, because it 
translates into a rather limited number of possible stable objects of any size. Ba- 
sically, there is nothing between atoms bound together by electromagnetism (and 
consisting almost entirely of nuclei bound by strong interactions) and Jupiter-sized 
brown dwarfs that can play a significant role in cosmology, the one exception be- 
ing primordial black holes which, above a mass of about 10 10 kg, would evaporate 
(Hawking 1974) sufficiently slowly to play an important dynamical role in the uni- 
verse. 

We can now proceed to describe some of the failures of a purely standard model 
early universe based solely on general relativity (SMGRC). 

Let us start with the implications of the relatively simple process of nucleosyn- 
thesis in the early universe alluded to above. This process, which takes place when 
the universe is just a few minutes old in the SMGRC (i.e. a few minutes after it 
"emerges" from the Big Bang singularity), results in the binding of about one quar- 
ter of all baryons into helium-4 nuclei, and the production of trace amounts of deu- 
terium, helium-3, and lithium-7. Relatively simple numerical calculations (Alpher 
1948; Peebles 1966ab, Wagoner 1967) show that the relative abundances of these 
isotopes are sensitive functions of the number density of baryons (compared to say 
the number density of photons). These isotopes can be identified spectroscopically 
in many objects in the contemporary universe - stars, HII regions, ... Although 
the current abundances have been changed to greater or lesser extents in different 
systems from the original primordial abundances, by chains of inference of vary- 
ing elaborateness and reliability, these primordial abundances can be inferred. A 
great success of SMGRC is that these several primordial abundances are mutually 
consistent, even though there is only one free parameter of much influence - the 
baryon-to-photon ratio. This turns out to be approximately 10~ 10 . 

Where success turns to disaster is when we combine the inferred baryon-to- 
photon ratio with the current photon density (the number density of cosmic mi- 
crowave background (CMB) photons, to be specific) which we can measure with 
sensitive microwave detectors (Penzias 1965). This allows us to infer the current 
number density of baryons, and hence the current density of Standard Model mat- 
ter, 

p SM * 5 x 10- 31 g/cm 3 . (2.1) 

The difficulty comes when this is compared to the total energy density when inferred 
from the rate of expansion of the universe using the Freedmann equation: 

h2= 8ttG M (22) 
3 

where H is the Hubble parameter - the logarithmic time derivative of the rate of 
the cosmic scale factor a(t). Since today H = H a ~ 70km/s/MpC, one finds that 

P = Pc = ^^lxlO- 29 g/cm 3 . (2 .3) 
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We thus see that fi& = Psm/Pc = 0.05. Within the context of SMGRC this discrep- 
ancy must mean that the geometry of the universe is strongly negatively curved 
since actually 

2 _ 87rGj\rp k 
H °-^^ + ^tf- (2 ' 4) 
Unfortunately, detailed observations of the power spectrum of fluctuations in the 
CMB (Lange 2001, Balbi 2000, Spergel 2003) set an upper limit on n curv = 1 - 
Psm/Pc < 0.05 which is significantly smaller than the inferred value in SMGRC of 
0.95. 

The only possibility within the SMGRC of resolving this missing energy-density 
problem is that the other 95% of the energy density of the universe is in the form 
of so-called primordial black holes (PBH) - in particular black holes formed before 
nucleosynthesis. There is no mechanism within the SMGRC for the formation of 
long-lived black holes. That requires a first order phase transition late enough after 
the singularity that the mass contained within the horizon can form a black hole 
with a lifetime to Hawking evaporation greater than the current age of the universe 
- about 10 billion years old. The only Standard Model phase transitions that qualify 
are the electroweak phase transition, the colour-confinement phase transition and 
the chiral symmetry breaking phase transition. But none of them are first order 
within the standard model. Furthermore, the universe can be shown (Carr 2010) 
not to contain black holes with twenty times the average mass density as the mat- 
ter, indeed Qpbh = Ppbh/Pc is constrained for Mpbh to be less than (at most) 
10~ 4 . Only for PBHs with mass less than or equal to the Planck mass is this limit 
inapplicable, but that is precisely because they are in the ultra-violet regime where 
GR meets quantum mechanics, and so the behaviour of such PBHs is ill determined 
in other words, they are outside the SMGRC. Furthermore, as we shall discuss 
below, having all of that 95% of the missing energy-density be of a form that scales 
with redshift (z) as (1 + z) 3 (as PBHs do) is not viable. Most of it needs to depend 
much more weakly on z. 

This was, of course, not the first case made for missing matter. The original 
inference that the matter that we see is not sufficient to source the gravitational 
forces that we infer from observations dates back to Zwicky (1933) and his observa- 
tions of clusters. Indeed, in all gravitationally bound systems from the length scale 
of dwarf galaxies up (ie. galaxies, groups and clusters) one observes that (except 
in the systems' central regions) the velocities of the objects that are apparently 
bound to the system exceed the local escape velocity in the gravitational potential 
that one would infer from the observed mass distribution using GR (i.e. Newtonian 
gravity). Thus, in spiral galaxies one observes that rotation curves are flat functions 
of radius outside the region where the luminous mass is concentrated, instead of 
falling as r -1 ' 2 as Newtonian gravity predicts. Similarly, clusters are observed to 
have X-ray gas that is far too hot to be bound in the Newtonian potential of the 
galaxies and intracluster gas that comprise the visible cluster. Indeed this missing- 
gravity problem has been verified in far too many systems for it to be appropriate 
to review here. It is perhaps noteworthy that this same problem does not apply to 
denser systems of the same mass as dwarf galaxies, such as globular clusters. 

This missing gravity problem also rears its head in the question of how the 
large scale structures that we currently observe in the universe formed. In the early 
universe, radiation dominates the energy density; however, the energy density of 
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radiation p r oc (1 + z) 4 and so matter eventually comes to dominate, even though 
the number density of photons is approximately 10 10 times that of baryons. During 
radiation domination, fluctuations in the energy density do not grow. Loosely, in 
the competition between growth of fluctuations under the influence of their own 
self gravity, and the damping of those fluctuations because of the expansion of the 
universe, the self-gravity loses. Once matter dominates, fluctuations do grow, but 
the density contrast increases only linearly as (1 + z). However, we have a snapshot 
of the universe - the CMB at rcdshift of 10 3 - and it shows fluctuations of only 
10~ 5 . It is thus impossible for non-linear structures such as galaxies and clusters to 
have grown in a SMGRC until today, much less by redshifts of a few when they are 
observed. Put another way - there is insufficient attractive self-gravity in SMGRC 
to explain the growth of the observed structures in the universe. 

The above failures of SMGRC are failures to provide sufficient ordinary at- 
tractive gravity, but attractive gravity alone is insufficient to explain the universe 
we see. In particular both the universe's current expansion and its expansion at 
some epoch in the very distant past seem to involve acceleration. The current 
accelerating expansion is inferred by direct observation of the recession velocity 
(rcdshift) of distant type la supernovae as a function of their distance, as derived 
from their apparent brightness. (Originally Riess 1998 and Perlmuttcr 1999. Most 
recently Amanullah 2010, Kowalski 2008, Kessler 2009, Contreras 2009, Balland 
2009, Bailey 2008 and Hicken 2009) (These supernovae are standard candles once 
one recalibrates to account for their different masses.) 

This accelerated expansion is also necessary to explain the age of our universe 
whose current hubble expansion rate is Hq — 72km/s/MpC. A flat matter domi- 
nated universe of this expansion rate is less than 10 billion years old, whereas the 
oldest globular clusters are thought to be 2-3 billions years older than that. An 
accelerating universe would be older than the objects in it. 

An earlier period of accelerated expansion is thought to be necessary to explain 
other conundrums about the universe in which we live. Most prominent among these 
is the homogeneity problem: if light of the CMB is just reaching here from after trav- 
eling through the universe since it was a brief fraction of its current age, then how 
is it that the CMB is a blackbody of temperature 2.7K in all directions? Shouldn't 
the emitting regions on opposite sides of the sky (say) be causally disconnected, 
and so at unrelated temperatures. This and a number of similar conundrums - the 
geometric flatness of the universe, its age, the high entropy - can all be addressed 
if the universe underwent a period of accelerated expansion, called inflation, early 
in its history. 

The attractive gravity among particles of matter or radiation acts to deccelerate 
the expansion of the universe. Within general relativity, acceleration is driven by 
forms of energy density that scales at least as slowly as (1 + z). The Standard 
Model has only one such form of energy density - the energy of the vacuum, p v . It 
is unclear within the Standard Model what is the value of p v . A naive estimate is 
approximately 10 120 times what is required to drive the current cosmic acceleration. 
This is known as the cosmological constant problem. Certainly the model does not 
contain energy density to drive both the early universe inflation, and the current 
acceleration. 

Finally, we mention briefly one other possible failure of SMGRC: the anomalous 
decceleration of the Pioneer 10 and 11 spacecraft. This anomalous acceleration 
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of ap = (8.74 ± 1.33) x 10~ 10 m/s 2 appeared in the orbit of the spacecraft once 
they entered the outer solar system and persisted until contact with them was 
left. It has long been claimed that no known force can provide the decceleration. 
Recently, a prosaic explanation has been suggested for at least part of the anomalous 
acceleration (Turyshev 2010), but that remains preliminary. 

3. The Standard Solution Three Dark Sectors 

The standard solution to this serial failure of the SMGRC is well known to the 
reader, and is known as the LambdaCDM (ACDM) concordance inflationary cos- 
mology. It involves introducing three entirely separate dark sectors: dark matter 
(non-relativistic particles), dark energy, and an inflaton (a second form of dark 
energy). 

The inflaton is a scalar field introduced to address the homogeneity and flatness 
problems through accelerated expansion in the early universe. This accelerated ex- 
pansion is the result of a period where the potential energy of the inflaton field 
is the dominant source of energy density and pressure in the universe. Because 
the pressure is negative, or equivalently because potential energy density does not 
dilute as the universe expands if the field is not evolving (no kinetic energy), the 
Friedmann equation 2.2 implies exponential expansion. The great success of the in- 
flationary model is that even as it explains the homogeneity and flatness, it provides 
a mechanism for generating fluctuations that later seed the growth of large scale 
structure. It is unfortunate that the natural amplitude of those fluctuations is 0(1), 
about five orders of magnitude larger than is observed. Taming those fluctuations is 
possible through mechanisms of varying (technical) naturalness that fine tune the 
parameters of the inflaton action. The disagreeable nature of this fine tuning has, 
over time, largely been forgotten given the success of the model in reproducing the 
Harrison- Zeldovich (scale free) spectrum of adiabatic perturbations that had previ- 
ously been hypothesized. The dominance of adiabatic fluctuations on a wide range 
of scales is consistent with inflation, though not a necessary feature of inflationary 
models. The phase coherence of those adiabatic modes evidenced by the observed 
doppler peaks in the CMB is characteristic of inflation. The moderately generic 
inflationary prediction that the spectrum deviates somewhat from scale free cither 
has been confirmed, or may soon be tested; similarly the generically predicted ten- 
sor perturbations are being sought, though their amplitude may be too small to be 
detectable. 

Efforts to connect the inflaton to Standard Model, or "just beyond the Standard 
Model" physics (eg. the Minimal Supersymmetric Standard Model) have not been 
successful, so that the range of candidate inflatons is enormous but with little or 
no foreseeable connection to physics that can be tested in the laboratory. 

The standard cosmology requires potential energy to dominate not just in the 
early universe, but again today. Today's dark energy density's properties are little 
better determined than those of the inflaton - it must be approximately three- 
quarters of the critical density or about (10 _3 eV) 4 - many many many orders of 
magnitude smaller than the inflaton energy density (except, see Dutta 2005). In the 
simplest model, it is just a cosmological constant, aka the vacuum energy density, 
with isotropic pressure, an equation of state w = p/p = — 1 and no anisotropic 
stress. All observations to date are consistent with this simplest model, although 
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that does not deter theorists from exploring "richer" (i.e. more complicated) models, 
nor observers from seeking to constrain those models. 

Finally (though historically first in priority), standard cosmology requires the 
introduction of dark matter, in particular non-relativistic (a.k.a. cold) dark matter. 
Numerous candidates have been identified. They range in mass from fractions of 
an electron volt (eg. axions, which are cold because they are produced as a bose 
condensate) to many solar masses. Some are ad hoc (i.e. invented exclusively or 
at least principally to solve the missing matter problem, with no connection to 
the standard model, nor to extensions of it that were devised to explain particle 
physics problems. Others are better motivated - for example the lightest supersym- 
metric partner (LSP) or the axion. The first great success of dark matter is that 
it addresses multiple problems simultaneously. With sufficient dark matter appro- 
priately distributed the rotation curves of spiral galaxies, velocity dispersions of 
elliptical galaxies and kinematics (including X-ray temperatures of clusters can be 
explained. The amount of dark matter that must be introduced is specific to each 
system, as is its precise distribution, but the variation is not unexpected. Beyond 
that, the existence of dark matter allows us to understand how large scale structure 
can grow by gravitational accretion even though the observed fluctuations in the 
CMB arc so small. This is because dark matter fluctuations start growing from the 
equality of the matter and radiation energy densities - earlier than baryonic matter 
fluctuations, which are limited by the coupling between matter and radiation un- 
til recombination. The dark matter fluctuations can therefore act as seeds for the 
growth of structure in the baryonic matter. 

The dark matter model made at least two further successes in two predictions 
that have since been confirmed. First, that light propagating past mass concentra- 
tions like clusters would be gravitationally lensed, and if one used GR to extract the 
mass (distribution) of the cluster, it would be consistent with the mass distribution 
as extracted by measurements of the kinematics of the constituent galaxies and 
gas (Tyson 1984). Second, that when structures dominated by dark matter (such 
as clusters) collided with each other, the dark matter halos (and stars) associated 
with each would scatter (essentially) only gravitationally, even as the gas associ- 
ated with each had more complex interactions. This prediction is consistent with 
observations of systems like the Bullet Cluster (Clowe 2006). 

We will keep these three predictions in mind as we explore the viability of 
alternatives to a multi-component dark sector. 

4. Modifying Gravity Instead 

The subtext of the previous section was that while our universe can apparently 
be modeled phcnomenologically by a multi-component dark sector, it cries out for 
a more parsimonious description. Adding dark matter or dark energy represents 
adding new forms of stress-energy T^ v , not included in the standard model, but 
preserving the Einstein field equations themselves - the dynamical equations in 
which matter sources gravity: 

G>„ = SttGnT^. (4.1) 

However, given that we have yet to find a model of the dark sector that elegantly 
wraps together dark matter, dark energy and inflaton, we might ask whether a mod- 
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ification of the dynamical equation itself (or more broadly of the underlying theory 
of GR) could give us a model of the Universe that is shaved more closely by Occam's 
Razor. Of course it would be advantageous if it was not only a simpler model, but 
could be experimentally or observationally differentiated from modifications of T M „ 
exclusively. 

There is at least some hint that such an effort at a unified theory of the devia- 
tions of the universe from SMGRC might succeed. That is that there is a common 
scale appearing in several of these deviations. The observed cosmic acceleration 
clearly sets a scale - c/Hq ~ 4000MpC. Alternately, this can be viewed as an accel- 
eration scale H n c = 7 x 10~ 8 cm/s . Interestingly, this is not so different from the 
typical value of the centripetal acceleration at which the rotation curves of galaxies 
deviate from the Newtonian (i.e. GR) prediction. Thus, the acceleration scale that 
MOdified Newtonian Dynamics (MOND), a phenomenological modified-gravity re- 
placement for dark matter (Milgrom 1983), introduces is an ~ 6H c. Amusingly, 
the anomalous acceleration of the Pioneer satellites (Anderson 2002) was also of 
this order, api oneer ~ 9 x 10 ?8 cm/s 2 . 

One can even perhaps squint and see hints of this scale in the inflationary 
theory. The inflationary prediction that all fluctuations in the universe are (seeded 
by) realizations of a Gaussian random statistically isotropic field, while consistent 
with observations at most angular scales, seem to fail at angular scales above 60° 
(Spergel 2003, Copi 2006). In particular the two point-correlation function of the 
CMB is close to vanishing on all angular scales between about 60 and 170 degrees. 
Because of cosmic variance, this is very unlikely if the coefficients of the spherical 
harmonics, ae m , are Gaussian random. Interestingly, 60° is the angle subtended by 
the length scale H^ 1 at z 1. 

The truth is that most modifications of gravity do not tie together all of these 
different failings of SMGRC, and so they are, at least in that sense, no better or 
worse than dark sector solutions. It may well be that instead of dark energy we have 
"accelerated expansion from gravity leaking to extra dimensions" (Deffayct 2001) as 
in DGP gravity (Dvali 2000), or a modest modification of the gravitational part of 
the Einstein-Hilbert action as in f(R) models (Carroll 2004 and references thereto). 

It is without a doubt a worthy goal to seek out modified gravity models of 
dark energy or inflation and look for points of observational differentiation from 
dark energy. But if the aim is to seek a model that provides a more comprehensive 
answers to the failures of the SMGRC, capitalizing perhaps on the coincidences of 
scale in the observed anomalies, then it will be necessary to address the replacement 
of dark matter by modified gravity. A preliminary framework for considering this 
possibility is MOND (Milgrom 1983). There is more than one version of MOND, 
some modifying Newton's inverse square law of gravitation (or Poisson's equation), 
some Newton's second law of motion (Milgrom 2005). We shall frame our discussion 
in terms of a modification of Poission's equation: 

V • V$|/ao)V$] = 4ttGp . (4.2) 

The function, fi(x) must be chosen to satisfy the necessary limits - namely that 
H(x) — > 1 for large x (the Newtonian limit) and fi{x) — > x as x — > (the MOND 
limit). Far from a static point mass, this MOND limit gives V$ — > s/GMao/r. For 
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circular orbits, this means that in MOND 

^circular -> (GMa ) 1/4 , (4.3) 

i.e. rotation curves are flat at large radii, as observed. Surprisingly, the form 

/i(x) = - 7 =L= (4.4) 

works very well for most galaxies. 

What does work well mean? First of all, 5.22 together with 4.4 reproduces the 
broad brush features of galactic rotation curves. Of course this is not surprising - 
it was chosen to give Newtonian gravity in the inner reaches where there are a lot 
of baryons, and flat rotation curves outside the baryons. But there is a little bit 
more to even this simple statement. After all, this this simple phenomenological 
model has one universal function ^i(x) and only one adjustable constant per galaxy 
the mass-to-light ratio. This is considerably less freedom than is available in 
fitting dark matter models to individual galaxies, as there one has an unknown 
dark matter density distribution (limited principally by expectations that such 
distributions be both reasonable (eg. non-singular) and reasonably similar among 
galaxies of similar morphology, and hopes that they could emerge from numerical 
modeling of galaxy formation. Moreover, the mass-to-light ratios that are derived 
for spiral galaxies from MOND are eminently reasonable - in other words they are 
around 20, consistent with expectations independent of MOND. 

Furthermore, MOND gives us a way to understand the well known Tully Fisher 
relation, that the luminosity of a spiral galaxy scales as the fourth power of its 
velocity width - the amplitude of its rotation velocity. After all 4.3 implies that 
M oc v 4 , so if mass-to-light ratios of spirals are reasonably uniform, then C oc v 4 
as well. Similar arguments (esp. Sanders 2010) have been made about the Faber- 
Jackson relation - the empirical relation between the luminosity of an elliptical 
galaxy and its velocity dispersion: C oc cr 4 . 

Even more surprisingly the agreement of MOND with galaxy rotation curves is 
not only at the qualitative level. Galaxy rotation curves have all sorts of bumps 
and wiggles. In a dark matter model, such features can always be attributed to 
features in the dark matter distribution. Not so in MOND - all features must be 
traced back to the baryon distribution, which alone sources the gravitational field. 
Amazingly (at least to the author), MOND succeeds in correctly reproducing the 
details of many (or perhaps all?) rotation curves from the details of the distribution 
of baryonic mass (Sanders and McGaugh 2002; McGaugh 2005). 

The situation for MOND at smaller scales is less clearly positive than at inter- 
mediate scales. Although dwarf galaxies and globular clusters are of similar mass, 
dwarf galaxies have a missing gravity problem, and globular clusters, it has been ar- 
gued, do not (Baumgardt 2009; Jordi 2009; Lane 2010; Conroy 2010). Of course, the 
mass of the system is not the only issue, especially when the system is not isolated, 
such as for most globular clusters embedded in the Galactic halo. The concentra- 
tion of the system also matters, as does the environment - especially at radii where 
the Newtonian acceleration falls below the MOND threshold. Thus others (Gentile 
2010) find that such globular clusters are fully consistent with MOND, and may 
even be demonstrably inconsistent with Newtonian gravity (McGaugh 2010). 
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At larger scales, it has long been known that (galaxy) clusters do not follow 
MOND. For example Aguirre et al. (2001) argued that "given observed gas density 
and enclosed mass profiles and the assumption of hydrostatic equilibrium, MOND 
predicts radial temperature profiles which disagree badly with observations." They 
showed this explicitly for the Virgo, Abell 2199 and Coma clusters, but argued 
that the results are general. The discrepancy could be resolved by the presence of 
additional non-luminous baryons, but it would have to be 1-3 times the observed 
cluster gas mass If this discrepancy is to be resolved by positing additional (presum- 
ably baryonic) dark matter, then this dark matter must at least equal the observed 
mass within 1 Mpc. Of course, one could just posit the presence of dark matter in 
the cluster. This may seem like a cop-out - the whole point of MOND is to avoid 
introducing new unknown forms of stress energy - but is not quite as bad as the 
first impression would suggest, because it could be hot dark matter. (Dark matter 
that was moving sufficiently fast at early times not to be captured into galaxies.) 
Why is this less distateful than cold dark matter? Because at least there are natural 
candidate - neutrinos. 

Although neutrinos are massless in the standard model, we know that neutrinos 
are not properly described by the standard model. In particular the eigenstates of 
flavour (i.e. electron, mu and tau neutrinos) are not the freely propagating eigen- 
states of the Hamiltonian. The canonical models for this neutrino mixing require 
that the neutrinos have non-vanishing masses. It had been suggested (Angus 2007) 
that a neutrino mass of 2 — 7eV would appropriately fill the need of clusters for 
dark matter in addition to MOND, however (Natarajan and Zhao) have since ar- 
gued against that model's sufficiency to explain observed weak lensing signals from 
clusters. Nevertheless, it continues to be argued (Feix) that, at least within mod- 
ified gravity theories with appropriate MOND limits, neutrino mass can play the 
role demanded. 

Finally, before we proceed to discuss some serious concerns, it should be noted 
that since the Bekenstein's introduction of TeVeS (Bekenstein 2004), there exist 
covariant theories based on an action principle that recover MOND (in its modified 
Poisson equation form) in the appropriate limits. TeVeS is a theory with. 

In Bekenstein's theory gravity is mediated by a tensor field g a b with associated 
connection W a b and inverse g ab , a one-form field A a , and a scalar field <j>. The 
Einstcin-Hilbert action governs the dynamics of g a b- 



(G is Newton's constant and R is the scalar curvature of g a b-) The dynamics of <j> 
are given by 



where we have used the convention of Skordis (2006). /i is a non-dynamical field 
and V is a free function that is chosen to give TeVeS a non-relativistic MONDian 
limit. The dynamics of A a are described by 
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where F a b — 2V[ a A^ and brackets denote antisymmetrization. Indices are raised 
with g ao and K is a dimensionless parameter. Variation with respect to the Lagrange 
multiplier held A enforces that A a be 'unit timelike': g ab A a At, = — 1. 

A key feature of TeVeS, aside from the breaking of Lorentz symmetry by the 
background time-like A a , is that matter couples to a different metric g a b than the 
dynamical one g ab described above. Test particles thus follow the geodesies of g not 
g. For some collection of matter fields f A the matter action is thus 



S m = J d 4 x(-g)iL[g,f A ,df 



When the unit-timelike constraint on A a is satisfied, the metrics, one-form field, 
and scalar field are related by the disformal transformation: 

g ab = e 20 (g ab + 2A a A b sinh(20)) .. (4.5) 

Zlosnik et al. showed (Zlosnik 2006) that TeVeS can be re-written as a Vector- 
Tensor theory akin to Einstein-Aether theories (see for example Eling 2004 for 
a review) with non-canonical kinetic terms. This enabled them (Zlosnik 2007) to 
broaden the class of covariant modified gravity models that have similar properties 
to TeVeS - appropriate non-relativistic MOND (and Newtonian limits). They all 
share the existence of a timelike vector field A a , a so-called Einstein Aether, and 
hence are referred to as Generalized Einstein Aether (GEA) theories. 

A subclass of these models are the so-called Generalized Einstein Aether GEA 
theories (Zlosnik 2007) that start with a covariant action and recover MOND (in 
its modified Poisson equation form) in the appropriate limits. In its simplest form, 
GEA contains a vector field, A coupled to Einstein-Hilbcrt gravity but not to 
matter: 



H +C(g,A) 



16ttG 



+ S M (4.6) 



where g a p is the metric, R the Ricci scalar of that metric. L is constructed to by 
generally covariant and local. For simplicity, the matter action, Sm, is taken to 
couple to g but not to A. 

As for TeVeS, A is enforced to be unit time-like through a Langrange multiplier. 
Otherwise, the Lagrangian is taken to depend only on covariant derivatives of A, 

C{A,g) = J^-T{K) + \{A a A a + l) 

IC af \ s = c l9 ^g ia + c 2 6^ + c 3 6^ (4.7) 

where q are dimensionless constants and M has the dimension of mass. (TeVeS is 
formally equivalent to a theory with an extended /C, and a more exotic method of 
achieving a non- vanishing background value for A.) More details can be found in 
(Zlosnik 2007). 

TeVeS and GEA possess important advantages over MOND alone. Most impor- 
tantly they are (at least nominally) comprehensive frameworks in which any ques- 
tion about classical gravity can in principle be answered. Thus it is appropriate to 
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ask questions about cosmology, stellar structure, compact objects, gravity waves, 
gravitational lensing, solar system dynamics. The answers may be difficult to obtain 
in practice, but unlike in MOND one cannot respond to disagreements with data by 
invoking the fact that the model strictly applies only to isolated/static/spherically 
symmetric systems. 

It is encouraging that within this covariant framework one can simultaneously 
realize both the non-relativistic MOND limit and a cosmology in which the scale 
factor evolution is consistent with observations - i.e. late-time acceleration, and 
inflation. (For GEA, see Zlosnik 2008; for TeVeS see Skordis 2005 and Diaz-Rivera 
2006.) One of the reasons that it is indeed relatively easy to accomodate both 
MOND on galactic scales and acceleration on cosmological scales is that they tend 
to probe different domains of the arbitrary functions that characterize the theory. 
For example, in GEA, for a universe described by the usual flat FRW metric, filled 
with a perfect fluid, 

K = 3 ^ ^ 

where, of course, H = ^, the dot denotes differentiation with respect to t, and, 
following a = C\ + 3c 2 + c 3 < 0. a is negative and hence so is JC. 
The modified 00 Einstein equation is (Zlosnick 2008) 

[l-J? K a]H* + \rM*=*!£-p (4.9) 

where Tk = dF /dJC. Restricting to T of the form 

J- = 7 (-/C) n . (4.10) 
equation 4.9 is more clearly a modified Friedmann equation: 



1 + e 



H \ 2(n-l) 

M 



H 2 = *¥p, (4.11) 



where 

e= (1 - 2n)7(-3a) n /6. 

We see that, for n = 1/2 the Friedmann equations are unchanged (e = 0). For 
n = 1 e = 7a/2 and the only effect (on scale factor evolution) is that Newton's 
constant is renormalized - G' = G/(l + e). For n = 0, we recover a cosmological 
constant, A ~ sign(—j)M 2 . More generally (Zlosnick 2008), there is late time 
accelerated expansion whenever 7(2n — 1) > 0. Moreover, not only can one obtain 
late time acceleration, but if n evolves as K evolves, then one could obtain a wide 
range of scale factor evolutions. 

A particular key point however is that cosmology probed only the K, < be- 
haviour of J ' . This is because it is sensitive to the time derivatives of A not its 
spatial derivatives. However, in the MOND limit, for example in a galaxy, the sys- 
tem is static. Thus it is is the spatial gradients of A a that dominate JC, and thus 
K > 0. 

It is in some ways unfortunate that structures and scale factor evolution probe 
disjoint ranges of JC. It means that it is difficult to promote such theories as offering 
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truly unified alternatives to dark matter plus dark energy. Rather one should say 
that both can be accommodated, and the fact that the relevant acceleration scale 
associated with each is roughly similar just means that it is easier to interpolate 
between the two domains of K. 



Having promoted the virtues of modifying gravity, it is time to get to the point of 
this talk - that it is hard to modify gravity and get away with it. We shall focus 
on a series of challenges that any attempt to replace the dark sector by modified 
gravity must face. 



One of the successes we outlined for dark matter is that it allows large scale 
structures to grow in the early universe, and without the dark matter, there is in- 
sufficient time to allow structure to grow. This is of course obvious, but there has 
always been the vaguely justified hope that, since MONDacts to strengthen gravity 
on large scales the increased attractiveness would result in more rapid structure 
growth. This of course ignores the role that stronger gravity plays in damping 
structure growth - faster growth necessitates that the attractive gravity work still 
harder to get structure to grow rather than be damped out by faster Hubble growth. 
Proving there is a no go theorem is not possible in the context of MOND because 
MOND is not well defined on cosmological scales. Nevertheless Lue and Starkman 
(Luc 2004b) were able to make some clear progress on this front relatively generi- 
cally. As they write: 

"By insisting that the gravitational interaction that accounts for the Newtonian 
force also drives cosmic expansion, one may kinematically identify which cosmolo- 
gies are compatible with MOND, without explicit reference to the underlying theory 
so long as the theory obeys Birkhoffs Theorem. Using this technique, we are able 
to self-consistently compute a number of quantities of cosmological interest. We 
find that the critical acceleration aO must have a slight source-mass dependence 
(aOM 1 / 3 ) and that MOND cosmologies are naturally compatible with observed 
late-time expansion history and the contemporary cosmic acceleration. However, 
cosmologies that can produce enough density perturbations to account for struc- 
ture formation are contrived and fine-tuned. Even then, they may be marginally 
ruled out by evidence of early [z < 20) reionization 

The calculation is reasonably straightforward. It generalizes the classic elemen- 
tary derivation of the Friedmann equation from conservation of energy. We repro- 
duce its outlines here. (For details sec Lue 2004a.) 

Consider a homogeneous universe filled with homogeneous matter of density 
p(t) <~ a~ 3 . Assume that it is described by the (flat FRW) line element 



with some specified scale- factor evolution, a(t). Gravitational dynamics must be 
altered in a specific way to achieve the observed cosmic expansion history, a(t). 
Writing a(t) as the solution to some alternative Friedmann equation: 



5. Troubles with Modifying Gravity 



(a) When Modified Gravity Means Dark Matter 



ds 2 = dt 2 - a 2 (t)5 lj dx l dx j , 



(5.1) 



a 2 /a 2 = H 2 g(x) 



(5.2) 
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(where x = ^wGp/H^ is a dimcnsionless density G is Newton's constant, and H is 
today's Hubble scale), we see that the function g(x) is determined by the observed 
a(t). If one requires that the fundamental gravitational theory respects Birkhoff's 
Theorem, then one can uniquely determine the metric of a spherically symmetric 
source (Lue 2004a). That metric is described by the line element 

ds 2 = g a(r)dt 2 - g rr {r)dr 2 - r 2 dO , (5.3) 

with g 00 (r) = g- 1 = 1 - r 2 H 2 g (r g /r 3 H 2 ) . (5.4) 

Here r g — 2GM, is the usual Schwarzschild radius of a matter source of mass M. 
Note that the metric components are completely determined by a(t). 

In MOND the gravitational acceleration exerted by a body of mass M obeys 

: -GMr- 2 \a\ > a 
\a\ < ao 

for some critical acceleration, a . To recover equation (5.5), we must have 

J x + cix 2 / 3 Einstein(x > x c ) 

9{X '^ \(3x 2 /Hnx + c 2 x 2 ^ MOND(Ki c ), ( ' 

for some constant parameters, (3, a, and c-i. Thus 

i-B |a|:-. 





where the critical MOND acceleration, ao, is 



9(3 2 (r g H ) 1 / 3 ] . (5.8) 



Observationally, we choose j3 ~ 15 so that for source masses the size of large galaxies 
(M ~ 10 11 M Q ), the critical acceleration is a ~ \H . To ensure that g(x) is 
continuous at a = a , C\ = c 2 + 3/3 [ln(3/3) — 1]. The remaining constant represents 
an arbitrary choice in zero-point energy for the Newtonian potential. Although 
d and c-i do not affect the Newtonian acceleration, they simulate curvature-type 
terms in the Friedmann equation despite the spatially-flat cosmology (see Eq. (5.1)). 
However, the resulting change in g rr (see Eq. (5.4)) causes only immeasurably small 
effects on gravitational lensing and the post-Newtonian parameter, 7. 

The form Eq. (5.7) for the MOND gravitational acceleration is slightly dif- 
ferent than that typically considered, where a = H /6 is a universal constant. 
The combination of homogeneous cosmology and Birkhoff's Theorem compels a 
weak dependence of a on the source mass, ao ~ M 1 / 3 . The confrontation of this 
prediction with data will be explored in (Kafka 2011). More immediately, this com- 
bination determines both the modified Friedmann equation, Eq. (5.2), and the full 
Schwarzschild-like metric of a spherical mass source, Eq. (5.4), without explicit 
reference to the details of the fundamental theory. With these two governing re- 
lationships, one may compute modifications of planetary ephemeris, gravitational 
lensing, growth of density perturbations , and determine the late-time integrated 
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Sachs- Wolfe (ISW) effect on the CMB. Details of these calculations for arbitrary 
g(x) are given in (Lue 2004). 

Here we recall particularly on accommodating late-time acceleration into MOND, 
and on the growth of fluctuations. Indeed, to obtain late time acceleration we must 
modify the form of g(x) at very low x: 

!x + 3/3x 2 / 3 [ln(3/3) - 1] x > (3/3) 3 (Einstein), 
/3x 2 / 3 ln(l + x) 0.1 < x < (3/3) 3 (MOND), (5.9) 

n A x<0.1 ("DarkEnergy"), 

where g(x) — £1^ w 0.7 is equivalent to a cosmological constant. Interestingly, as 
discussed in (Lue and Starkman 2004), if (3 were an order-of- magnitude larger or 
smaller, one could not extend the MOND regime all the way to the deceleration- 
acceleration transition and still be able to maintain both H ~ Hq as well as a/a ~ 
Hq. On the other hand, this is a rather complicated form for g(x), and there is no 
particular reason to expect this to emerge naturally from a covariantized MOND 
theory. 

Equations (5.2) and (5.9) allow us to compute the growth of linear perturbations 
in this cosmology. In the absence of dark matter, we must wonder whether MOND 
has sufficient growth from recombination to produced the observed structures at 
z ~ few. As discussed in (Lue and Starkman 2004), the evolution of linear density 
perturbations can be obtained in closed form. Consider a uniform overdensity in a 
localized spherical region such that 

p(t) = p(t) [1 + S(t)} , (5.10) 

where p is the background matter density that follows cosmological evolution. Pa- 
rameterizing time-evolution using x — 8ttGp/3Hq, the growing perturbation mode 
5(x) goes as (Lue 2004 and Multamaki 2003): 

'w^wfiwr ,5 - u) 

The normalization is fixed by requiring that 5 <~ few x 10~ 5 at recombination. In 
pure matter-domination for Einstein FRW, 5 oc x^ 1 / 3 oc a(t). Even if growth were 
as large as this the growth of density perturbations would be insufficient to account 
for the observed structure formation. But in each of the three regimes in Eq. (5.9), 
growth is slower than that given benchmark, 5 < a(t). It may seem counterintuitive 
that growth is suppressed during the MOND regime, when the self-gravitation of 
overdensities is enhanced, but the stronger gravity drives faster expansion. This 
suppresses perturbation growth more than the stronger self-gravity enhances it. 

This suppression of structure growth would seem to pose a significant difficulty 
for cosmological incarnations of MOND. One might point to the fact that in "Real 
MOND" ao does not scale with M, and so hope that these problems are peculiar to 
the above calculation. However, since that relied only on the validity of Birkhoff 's 
Theorem, which we shall argue below one wants to relinquish as gently as possible 
if at all, it seems unlikely that this conclusion can be avoided without significant 
consequences. 
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In light of the previous argument, it was surprising (at least to the author) to 
learn that in TeVeS and GEA the growth of structure is not suppressed. However, 
the very way that covariant versions of MOND evade this constraint is enlightening 
and arguably undermines the MOND program. Skordis et al. (2003) were the first 
to study the growth of perturbations in TeVeS. Their numerical analysis identified 
growing modes in the TeVeS scalar field, which support the growth of fluctuations 
in the Newtonian potential. Dodelson and Liguori (2006) showed analytically that 
in fact it was growing modes in the vector fields that were responsible for the growth 
of structure in TeVeS. Zlosnik et al. (2004) showed the same thing for GEA - a 
vector instability could replace dark matter as a seed for the growth of baryonic 
structure after recombination. 

We thus seem driven to the conclusion that MOND - which was designed to 
avoid the addition of a new undiscovered degree of freedom that source the gravi- 
tational field in galaxies and clusters today and in the linear perturbations that are 
their precursors,, requires instead the introduction of a new undiscovered degree 
of freedom that source the gravitational field in those linear-perturbations. This 
is in addition to the unexpectedly large neutrino masses that must be introduced 
to explain the dynamics of clusters. In the light of results of Lue et al. (2003 and 
2004), it is difficult to see how this problem can be made to away. 

Much has been made of the failure of MOND to model the Bullet Cluster (as well 
as other non-static, non spherically symmetric systems), which has been taken as "... 
direct empirical proof of the existence of dark matter" (Clowe 2006 ). Specifically, in 
these systems the gravitational weak lensing centers are coincident with the stars, 
but not with the hot x-ray gas. This is what one would expect from a dark matter 
model. If two clusters collide, the dark matter halos, and the stellar distributions, 
should pass through each other, perturbing each others' paths only slightly, through 
gravity. Meantime, the hot x-ray gas from the two colliding clusters appear to have 
interacted non-gravitationally, and therefore slowed down relative to the stars. 

More accurately, such systems might be portrayed as confirming the predictions 
of the cold dark matter theory. However, they should not have been unexpected 
in the context of MOND. After all, as we described, clusters are known, even in 
MOND, to require some additional source of gravity both today and at the time 
of their formation. Why then would we expect the weak lensing signal to coincide 
with the mass distribution? Thus it is probably reasonable that one can find ways 
to recover the observed weak lensing signal of the bullet cluster in the context of 
modified gravity theories (Dai 2008; Brownstein 2007). As Dai writes regarding 
GEA "As vector-field fluctuations are responsible in GEA for seeding baryonic 
structure formation, vector-field concentrations can exist independently of baryonic 
matter. Such concentrations would not be expected to be tied to baryonic matter 
except gravitationally, and so, like dark matter halos, would become separated from 
baryonic matter in interacting systems such as the bullet cluster. These vector- 
field concentrations cause metric deviations that affect weak lensing. Therefore, the 
distribution of weak lensing deviates from that which would be inferred from the 
luminous mass distribution, in a way that numerical calculations demonstrate can 
be consistent with observations. This suggests that MOND-inspired theories can 
reproduce weak lensing observations, but makes clear the price: the existence of 
a coherent large-scale fluctuation of a field(s) weakly tied to the baryonic matter, 
not completely dissimilar to a dark matter halo." On the other hand, it is equally 



Article submitted to Royal Society 



Modifying Gravity: You Can't Always Get What You Want 17 



plausible that a specific modified gravity theory would require the addition of cluster 
dark matter, as has been argued for TeVeS (Fcix 2007). 

But whether the missing gravity of clusters is caused by hot dark matter, or 
some dark field halo, it is clear that theories with MOND limits need to have physics 
beyond MOND - dark degrees of freedom - to explain clusters. 



(b) Keeping cosmology out of your galaxy 

One difficulty for covariant versions of MOND that has not yet been well ex- 
plored is the tension between the MOND limit and cosmology (but see Jacobs 2011) 
Recall from above that the (or at least a possible) starting point of MOND is to 
write a modified Poisson equation: 

^ (|V$|/a ) V$ =AttG P . (5.12) 

At least this is meant to be the correct non-relativistic static spherically symmetric 
limit of the correct underlying covariant theory. It is meant to apply to galaxies 
and clusters. But cosmic structures exist within a cosmic background which is de- 
cidedly not static. Thinking of MOND naively in terms of (5.12), one is inclined 
to ignore cosmology, as the "Hubble acceleration" is roughly Hfir, where H is the 
Hubble constant today and r is the distance to the center of the gravitating object 
(e.g. galaxy or cluster). The same cannot be said, in general, for a corresponding 
covariant theory. In GEA for example (the same is also true of BIMOND) invariant 
scalars are built out of derivatives on the metric(s) and, in particular, this means 
that time derivatives are put on the same footing as gradients. Due to cosmolog- 
ical considerations, therefore, one will need to compare the natural value of time 
derivatives (~ Ho) with spatial derivatives of gravitational potential(s). Since the 
MOND acceleration scale ao is of order H , serious tension is possible. 

For example, consider two different metrics to approximate that of a galaxy, one 
static: 

5oo = -(1 + 2$) (5.13) 
Qij = (1 - 2*)ay (5.14) 

and the other dynamical (accounting for cosmological expansion) : 

5oo = -(1 + 2$) (5.15) 
9ij = (1 - 2$>(t) 2 % (5.16) 

where, in both metrics, we assume goi = 0. If we choose the static metric, (5.13) 
and (5.14), we find the Poisson equation in GEA to be modified as 

V • [(2 + ciJ 7 ') V$j = 8irGp (5.17) 

with the scalar JC given by 

K=^pt + G { &) (5.18) 
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Thus, in order to recover (5.12) in the limit V$ <C M, one finds the free function 
to take the form 

T(JC) ~ aK + P>C 3/2 (5.19) 

for some constants a and (3. 

Considering we do live in an expanding universe, however, it is arguably more 
correct to use a perturbed FRW metric. While our experience tells us to expect that 
cosmology should be irrelevant when investigating e.g. solar system dynamics, at 
large separations it is clearly present, therefore nothing should stop us from using a 
metric which accommodates both limits. However, using the dynamic metric, (5.15) 
and (5.16), one finds (Jacobs 2011) that the modified Poisson equation is further 
modified: 

1 



M'f I ^7 • '!' ) •••!// '.:! ouF K ) 



+ V- (2 + Cl J)c)V$ = SwG{p + Sp) 



Also the scalar JC is now given by 
1 



K 



M 2 



3aH 2 (1 - 2$ + $ 2 ) - Cl I — 



(5.20) 



(5.21) 



But in the regime of greatest interest V<I> is of order H = H . Thus, it appears that 
Hubble expansion is relevant at the very point where the MOND regime begins. It 
is not clear, and in fact seems unlikely, that the F{K) given by (5.19) would still 
be sufficient for a successful theory of MOND. 

We now see what the unexpected price of "covariant-izing" MOND can be - 
the same mechanism that allows one to modify the Poisson equation, namely using 
tensors built from derivatives of the metric, puts Hubble expansion on the same 
footing as gradients of the Newtonian potential(s). Factors of H that are introduced 
are non- negligible since MOND is necessarily concerned with the limit where V$ is 
of order H . Cosmology apparently becomes inescapable luggage in this endeavor, 
an issue that would appear to be fairly generic in any attempt to create a covariant 
theory of MOND. Interestingly, Milgrom's BIMOND (Milgrom 2009) docs seem to 
be an exception. 



(c) Living without Birkhoff? 

The final point of concern that will be raised here is that MOND, and covariant 
theories with MOND as a limit generically lack Birkhoff 't theorem. 

In Newtonian gravity, the inverse square law guarantees that the gravitational 
flux through an enclosed surface is conserved, this is Gausss law. In General Rel- 
ativity Birkhoffs Theorem (BT) plays the same role. Two particularly important 
aspects of BT are that: (a) outside a spherically symmetric mass distribution, the 
metric depends on the distribution of the matter density only at second order in 
the potential, due to the role of binding energy as a source of gravity (and that 
serves only to change the mass of the system) ; and (b) a spherically symmetric shell 
of mass has no effect on the metric in its interior. We depend on these properties 
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to allow us to calculate almost any gravitational fields without knowing details of 
the distribution of matter all over the Universe. Indeed many calculations in grav- 
ity would become undoable either in principle or practice, and many others would 
become enormously more difficult. 

In general, modifications to GR tend to violate BT (Satz 2005; Moffat 2007). 
And yet, almost all calculations of the implications of those modifications rely on 
BT, or on other congenial averaging properties of GR that seem unexpected in a 
non-linear theory, and that one should not therefore assume persist in modifica- 
tions of GR. After all, even in GR the averaging problem remains unresolved - to 
what extent is the average of the Einstein tensor in an inhomogeneous universe the 
solution of the Einstein equations with the spatially averaged stress-energy as a 
source? 

It is important to appreciate that BT, or more specifically, the fact that a 
spherically symmetric shell exerts no gravitational field on its interior is a rather 
surprising feature of GR. In Newtonian physics, we can understand that for a mass 
at a point in the interior of the shell the pull from a spherical cap on one side is 
perfectly balanced by the pull from the spherical cap of the same angular size on 
the opposite side. This is true even for an off-center location. The same is true only 
for a test mass in GR. A real off-center mass distorts the geometry so that the 
shell is no longer spherical. This leads to a force on the mass toward the center! 
However the magnitude of the force is suppressed both by the ratio of the point 
mass to the shell mass, and by the ratio of the Schwarzschild radius of the shell to 
its actual radius. Thus, for most cosmologically relevant situations, the corrections 
to the Newtonian intuition are small. 

Now MOND too has a Birkhoff-like theorem. In the interior of a spherical shell, 

V • [M(|V$|/a )V$] = -S- V$ = 0! (5.22) 

The problem is that as soon as an off-center mass is introduced into the interior 
of the shell, the theorem does not apply, and, unlike in GR, the corrections are no 
longer small (Dai 2008). The reason is that it was not actually BT that protected 
us in GR, but rather the inverse-square-law. Because the whole point of MOND 
is that gravity falls off more slowly than r~ 2 , breaking the spherical symmetry 
immediately leads to sizable forces in the interior of a shell. 

Cosmologically and astrophysically, above the scale of a single star (and even 
then), we almost never precisely satisfy the requirement of spherical symmetry. 
Galaxies are not spherical. They arc embedded in groups or clusters that are not 
spherically symmetric. Those clusters are themselves in environments that break 
spherical symmetry. In the context of perturbation theory, the question may be 
phrased in terms of unexpected mode-mode coupling. 

The question then is how large are the corrections to the conclusions that we 
draw through the false application of our inverse-square-law-based intuition to sys- 
tems with a more long-range force. The answer to that question is still not clear. 
Some progress has been made on simplified systems - sets of spherical shells, shells 
with off-center point masses, .... For example (Dai 2008), we know that interior 
point masses that are O(10%) of the mass of the surrounding shell can experience 
"anomalous" forces of a few percent of the MOND acceleration. On the other hand, 
it has also been shown (Matsuo 2009) that a spherical shell screens the effects of a 
perturbation on a shell exterior to it, reducing the field to first order in the pertur- 
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bation (0(e)) to 0(e 1 ^ 2 ). More work clearly remains to be done to provide confident 
answers for realistic situations to questions like - can we smooth the mass density 
distribution of a collection of stars and gas in order to calculate the local field at 
some point interior to that distribution, and if so, what is the level of accuracy 
that that approximation affords? How does the answer change if it is a collection of 
galaxies and gas in a cluster? or clusters in the universe? Living with out the 1/r 2 
force law is a dangerous thing to do, especially in a bumpy universe. The universe, 
or it subsystems, may simply not be accurately computable. 

6. Summary 

A model based solely on General Relativity and the Standard Model of particle 
physics fails to account for numerous phenomena on scales from the size of our 
solar system up to the size of the observed universe. These include the kinematics 
of galaxies and clusters and the thermodynamics of their gas, the lensing of light 
by these structures, and their very existence so soon after the big bang. On smaller 
scales, the Pioneer satellites appear to have an anomalous sunward component to 
their acceleration (although this may or may not have finally been attributed to 
anisotropic thermal radiation from the satellites). On larger scales, the universe 
is unreasonably homogeneous, isotropic and flat, unexpectedly old given its rate 
of expansion, expanding surprisingly rapidly given the low abundance of matter, 
expanding at an accelerating rate, and distinctly lacking in angular correlations on 
the largest angular scales. The dimensionful scales (lengths, times, accelerations) at 
which of many of these anomalies appear seem to be similar: H times appropriate 
powers of c. 

The concordance model of the universe addresses many of these puzzles through 
the introduction of three separate dark sectors - dark matter, dark energy and the 
inflaton. It is a very successful model (though mysteries remain), but its arbitrari- 
ness motivate us to consider the possibility that it is not, or not just, that the 
stress energy content of the Standard Model is incomplete but that our theory of 
gravity needs modification. The most challenging sector of the concordance model 
to replace with modified gravity is dark matter, because it ties together with rea- 
sonable success many of the problems. Modified Newtonian Dynamics provides a 
phenomenological framework within which to address many of these questions. It 
has particular successes of its own, such as explaining the relationship between 
galaxy luminosities and kinematics. However it also forces us to give up rather 
useful features of GR like Birkhoff's theorem. 

Covariant theories have been built that approach MOND in the appropriate 
limits. These theories may address some of the challenges of MOND - such as the 
failure to explain clusters without the introduction of hot dark matter, and the 
failure to grow large scale structures. However they have their own challenges. 

The program of exploring modified gravity continues to be a worthwhile one, 
and most excitedly modifying gravity on cosmological scales to explain cosmic ac- 
celeration seems both doable and testable. Replacing dark matter with modified 
gravity seems a more difficult road to travel despite the construction of covariant 
theories with MOND limits. At least we discover that they only seem to work if they 
bring back dark matter in the form of dark fields. That may be how the Universe 
works, but it isn't what the doctor (Milgrom) ordered. 
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